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Whistler Wave Cascades in Solar Wind Plasma 
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ABSTRACT 

Nonlinear three dimensional, time dependent, fluid simulations of whistler wave tur- 
bulence are performed to investigate role of whistler waves in solar wind plasma tur- 
bulence in which characteristic turbulent fluctuations are characterized typically by 
the frequency and length scales that are respectively bigger than ion gyro frequency 
and smaller than ion gyro radius. The electron inertial length is an intrinsic length 
scale in whistler wave turbulence that distinguishably divides the high frequency solar 
wind turbulent spectra into scales smaller and bigger than the electron inertial length. 
Our simulations find that the dispersive whistler modes evolve entirely differently in 
the two regimes. While the dispersive whistler wave effects are stronger in the large 
scale regime, they do not influence the spectral cascades which are describable by a 
Kolmogorov-like spectrum. By contrast, the small scale turbulent fluctuations 

exhibit a Navier-Stokes like evolution where characteristic turbulent eddies exhibit a 
typical k~^/'^ hydrodynamic turbulent spectrum. By virtue of equipartition between 
the wave velocity and magnetic fields, we quantify the role of whistler waves in the 
solar wind plasma fluctuations. 

Key words: (magnetohydrodynamics) MHD, (Sun:) solar wind, Sun: magnetic fields, 
ISM: magnetic fields 



1 INTRODUCTION 

The solar wind is an excellent in-situ laboratory for inves- 
tigating nonlinear and turbulent processes in a magnetized 
plasma fluid since it comprises a multitude of spatial 
and temporal length-scales associated with an admixture 
of waves, fluctuations, structures and nonlinear turbu- 
lent interactions. The in-situ spacecraft measurements 
(Matthaeus & Brown 1988, Goldstein et al 1995) reveal 
that the solar wind fluctuations, extending over several 
orders of magnitude in frequency and wavenumber, describe 
the power spectral density (PSD) spectrum that can be 
divided into three distinct regions (Goldstein et al 1995, 
Leamon et al 1999). The frequencies, for instance, smaller 
than 10-5 Hz lead to a PSD that has a spectral slope of 
-1 . This follows the region that extends from 10-5 Hz to 
or less than ion/proton gyrofrequency where the spectral 
slope exhibits an index of -3/2 or -5/3. The latter, a 
somewhat controversial issue, is characterized essentially 
by fully developed turbulence and can be followed from 
the usual magnetohydrodynamic (MHD) description. The 
spacecraft observations have further revealed that length 
scales beyond the MHD regime, that are smaller than 
ion gyro radius {kpi S> 1) and temporal scales bigger 
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than ion cyclotron frequency u > ujd = eBo/m^c, (where 
k, pi,ujci,e, Bo,me, c are respectively characteristic mode, 
ion gyroradius, ion cyclotron frequency, electronic charge, 
mean magnetic field, mass of electron and speed of light) 
exhibit a spectral break where the inertial range slope of the 
solar wind turbulent fluctuations varies between -2 and -5 
(Smith et al 1990, Goldstein et al 1994, Leamon et al 1999). 
Notably, the dynamics of the length-scales in this region 
cannot be described by the usual MHD models that possess 
characteristic frequencies smaller than ion gyro frequencies. 
At 1 AU, the ion inertial leng th scales are smaller than 
ion gyro radii in the solar wind (|Goldstein et al 19951 ). The 
latter, associated with plasma motion due to finite Larmor 
radii, can readily be resolved in the usual MHD models by 
introducing Hall terms to accommodate the ion gyro scales 
up to scales as small as ion inertial length scales. The higher 
time resolution databases identifying the spectral break 
indicate that Alfvenic MHD cascades (Smith et al 1990, 
Goldstein et al 1994, Leamon et al 1999) are terminated 
near the spectral break. The characteristic modes in this re- 
gion are observed to evolve typically on timescales involving 
the dispersive kinetic Alfvenic fluctuations. The onset of the 
second or the kinetic Alfven inertial range is still elusive to 
our understanding of the solar wind turbulence and many 
other nonlinear interactions. Specifically, the mechanism 
leading to the spectral break has been thought to be either 
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mediated by the kinetic Alfven waves (KAWs) (Hasegawa 
1976), or by electromaRnetic ion-cyclotron- A lfven (EMICA) 
waves taarv et al 20081 : IWu fc Y oon 2007*1. or by a class 
of fluctuations that can be dealt within the framework 
of the HMHD plasma model (Alexandrova et al 2007, 
2008; Shaikh & Shukla 2008, 2008a). Stawicki et al (2005) 
argue that Alfven fluctuations are suppressed by proton 
cyclotron damping at intermediate wavenumbers so the 
observed power spectra are likely to consist of weakly 
damped magnetosonic and/or whistler waves which are 
dispersive unlike Alfven waves. Moreover, turbulent fluc- 
tuations corresponding to the high frequency and kpi 2> 1 
regime lead to a decoupling of electron motion from that 
of ion such that the latter becomes unmagnetized and 
can be treated as an immobile neutralizing background 
fluid. While whistler waves typically survive in the higher 
frequency (and the corresponding smaller length scales) 
part of the solar wind plasma spectrum, their role in 
influencing the i nertial range turbu l ent spectral cascades 
is stiU debated llBiskamp et al 199^: IShaikh fc Zank 200S : 



Shaikh & Zank 2005; IShaikh fc Shukla 2005 

Shaikh & Shukla 2008). The Kolmogorov like dimen- 
sional arguments indicate that propagation of whistlers in 
the presence of a mean or an external constant magnetic 
field may change the spectral index of the inertial range tur - 
bulent fluctuations from to (Biskamp et al 19961). 

By contrast, the num e rical simulations (Biskamp et al 199Q ; 
IShaikh fc Zank 20031 : IShaikh fc Zank 20051 ) suggest that 
whistler waves do not influence the spectral migration of 
turbulent energy in the inertial range despite strong wave 
activity and that the turbulent spectra corresponding to 
the electron fluid fluctuations in whistler wave turbulence 
continue to exhibit a Kolmogorov-like fc"^''"' spectrum. 
What is not clear from these work is the qualitative role 
of whistler and the corresponding mode coupling interac- 
tions that mediate the inertial range turbulent spectra. 
Furthermore, t he whistler wave turbulence described 
in the Refs. llBiskamp et al 199^ : IShaikh fc Zank 20031 : 
IShaikh fc Zank 20051 ) focus on purely two dimensional 
interactions, ignoring thus the variations in the third 
dimension. It is therefore unclear whether the nonlinear 
whistler mode coupling interactions in three dimensions 
modify the energy cascades in the inertial range turbulence. 
Intrigued largely by these issues, the primary goal of this 
paper is to investigate the nonlinear interaction amongst 
whistler waves and turbulent fluctuations, based on nonlin- 
ear fluid simulations, in w > Ud regime where correlation 
length scales of turbulence are comparable to the electron 
inertial length scales. Understanding the role of whistler 
waves in the turbulent cascades is crucial to many other 
space plasma processes since whistler waves, in addition 
to solar wind turbulence, are instrumental in governing 
nonlinear processes in numerous other plasma systems 
that range fro m solar wind pCrafft fc Volokitin 2003 ; 



Saito et al 2"ooi : IStawicki et al 20051 : iGarv et al 2~ 
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Bhattachariee et al 199^ ). magnetic reco nnection 



the Earth's magnetosphere l|Wei et al 20071 ) to interstel- 
lar medium (|Burman 19751 ) and astrophysical plasmas 
l|Roth 20071 ) where characteristic fluctuations can typically 
be of several astronomical units. These are only a few 
of the numerous other studies. For more literature, the 



read ers can refer t o the simulation work by Biskamp 
(Bis kamp et al 19961 ) and others including Shukla (1978), 
Shukla et al (2001), Cho & Lazarian (2004), Gahier (2008), 
Urrutia et al (2008), Saito et al (2008), Bengt fc Shukla 
(2008), Shaikh (2009) and numerous references therein. 

In this paper, I focus on understanding the nonlinear 
turbulent cascades mediated by whistler waves in a fully 
three dimensional geometry. Our objective is to investigate 
the role of whistlers in establishing the turbulent equipar- 
tition amongst the modes that are responsible for the non- 
linear mode coupling interactions which critically determine 
the inertial range power spectra. Remarkably, we find that 
despite the equipartition processes mediated by whistler 
modes for which the wave activity is strong, the inerital 
range spectra continue to exhibit a Kolmogorov-like spec- 
trum where whistler wave effects are unimportant. We begin 
in Section 2 by describing the underlying whistler wave tur- 
bulence model and it's linear properties. Section 3 describes 
nonlinear simulation results of inertial range turbulent spec- 
tra. In section 4, we discuss the theoretical arguments cor- 
responding to the whistler wave turbulent spectra that cor- 
respond to the characteristic length scales smaller as well 
as bigger than electron inertial length (de). The process of 
turbulent equipartition between the magnetic and velocity 
field, quantifying the whistler wave effects, is also desribed 
in this section. Finally, section 5 summarizes our results. 



2 WHISTLER WAVE TURBULENCE MODEL 

Whistler modes are excited in the solar wind plasma when 
the characteristic plasma fluctuations propagate along a 
mean or background magnetic field with frequency lj > Uci 
and the length scales are c/u}pi < £ < c/ujpe, where ujpi,ujpc, 
are the plasma ion and electron frequencies. The electron 
dynamics plays a critical role in determining the nonlinear 
interactions while the ions merely provide a stationary 
neutralizing background against fast moving electrons and 
behave as scattering centers. The whistler wave turbulence 
can be described by the e lectron magnetohyd rodynamics 
(EMHD) model of plasma (i Kingsep et al. 19901 ) that deals 
with the single fluid description of quasi neutral plasma. The 
EMHD model h as been discussed in c onsiderable detail in 
earlier work llKingsep et al. 1990l : iBiskamp et al 199a : 
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Dastgeer et al. 2000al : iDastgeer et al. 2000b ; 

Shaikh fc Zank 20031 : IShaikh fc Zank 20051 ). In whistler 
modes, the currents carried by the electron fluid are im- 
portant, and we therefore write down only those equations 
which are pertinent to electron motion. These are electron 
fluid momentum, electric fleld, currents, and electron 
continuity equations, 

men— -^-l-Vs-VVe = -enE xB-VP-/imenVe, (1) 

ot c 
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The remaining equations are B = V x A, J = — enVg, V ■ 
B = 0. Here me,n,Ve are the electron mass, density and 
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fluid velocity respectively. E, B respectively represent elec- 
tric and magnetic fields and 0, A are electrostatic and elec- 
tromagnetic potentials. The remaining variables and con- 
stants are, the pressure P, the collisional dissipation /x, the 
current due to electrons flow J, and the velocity of light 
c. The displacement current in Ampere's law Eq. Q is ig- 
nored, and the density is considered as constant throughout 
the analysis. The electron continuity equation can therefore 
be represented by a divergence-less electron fluid velocity 
V • Ve — 0. The electron fluid velocity can then be associ- 
ated with the rotational magnetic field through 

-V X B. 



(5) 

On taking the curl of Eq. ([1} and, after slight rearrangement 
of the terms, we obtain 



^ - Ve X (V X P) + 
where 



(6) 
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and ^ = -meVe 
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Here P is generalized electron momenta. The curl of Eq. 
([6} eliminates the gradient of the scalar quantity (the third 
term from the left in the Ihs) and yields 



— + V X (Ve X f2) = -^me V X Ve 

where 



(7) 



n = V X P 



It can be seen from Eq. (6) that in the ideal whistler mode 
turbulence (i.e. neglecting the term associated with the 
damping /i), the Curl of generalized electron momenta is 
frozen in the electron fiuid velocity. This feature is strik- 
ingly similar to Alfvenic turbulence whe re the magnetic field 
is frozen in the ideal two fluid plasma l|Biskamp 20031 ). On 
substituting $7 into the above equation and using appropri- 
ate vector identities, we obtain the three-dimensional equa- 
tion of EMHD describing the evolution of the magnetic fleld 
fluctuations in whistler wave, 



d_ 

dt 



(B - d^V^B) + Ve ■ V(B - dlV^B) - 
(B - deV^B) ■ We = ^deV^B. 



(8) 



The length scales in Eq. ([8)l are normalized by the electron 
skin depth de = c/upe i.e. the electron inertial length scale, 
the magnetic field by a typical amplitude Bo, and time by 
the corresponding electron gyro-frequency. In Eq. the 
diffusion operator on the right hand side is raised to 2n. Here 
n is an integer and can take n = 1, 2, 3, • • •. The case n = 1 
stands for normal diffusion, while n = 2, 3, ■ ■ ■ corresponds 
to hyper- and other higher order diffusion terms. 

The linearization of Eq. ([8]) about a constant magnetic 
field B = Boz+B, where Bq and B are respectively constant 
and wave magnetic fields, yields the following equation. 



ajfc(l + d^fe^)B + 



CBq 



ifciik X B = 0. 



(9) 
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On eliminating the wave perturbed magnetic field from the 
above relation, one obtains the following dispersion relation. 



(10) 



where uJcq ~ eBo/mc and fc^ = k'^ + ky. The use of Eq. (|10|l 
in Eq. (9) leads to the following relation between the wave 
magnetic field and the velocity field. 
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±-k X B 
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(11) 



The rhs of Eq. (11), in combination with Eq. ((5]), corre- 
sponds essentially to the whister wave perturbed velocity 
field. This equation indicates that whistler waves are trans- 
verse and are produced by rotational magnetic field that 
leads essentially to the velocity field fiuctuations. On re- 
placing the rhs in Eq. (11) with the perturbed velocity field, 
it can be shown that the whistler modes obey equiparti- 
tion between the magnetic and velocity field components as 
fc^lBp ~ |Vep. The whistler wave activity can thus be quan- 
tified by how closely the characteristic modes obey the tur- 
bulent equipartition relation. In section 4, we will investigate 
using nonlinear 3D fluid simulations the equipartition medi- 
ated by whistler waves in the inertial range turbulent spectra 
to signify the role of whistlers in the solar wind plasma. 

It becomes evident from Eq. ^ that there exists an 
intrinsic length scale de, the electron inertial skin depth, 
which divides the entire turbulent spectrum into two re- 
gions; namely short scale (fcde > 1) and long scale (fcde < 1) 
regimes. In the regime fcde < 1, the linear frequency of 
whistlers is cok ~ kyk and the waves are dispersive. Con- 
versely, dispersion is weak in the other regime fcde > 1 since 
Wfc ~ fcj,/fc and hence the whistler wave packets interact 
more like the eddies of hydrodynamical fluids. The equation 
of EMHD [Eq. (|8])] is also exactly integrable, yielding the 
total energy integral; 



£=2/(1 



-d V^)|B|^dV 



and generalized helicity. Here d v is a 3D volume element. 
In the presence of dissipation (fi) the total energy decays 
eventually with time since 



dE 2 



d\ [{V'^Bf 



(12) 



Hence the inclusion of dissipation will damp the smaller scale 
fluctuations in the whistler wave turbulence. The damping 
of the smaller dissipative scales are not expected to influence 
the inertial range turbulent cascades. 



3 SIMULATION RESULTS 

Turbulent interactions mediated by the coupling of whistler 
waves and inertial range fluctuations are studied in three 
dimensions (3D) based on a nonlinear 3D whistler wave 
turbulence code that we have developed at Center for 
Space Plasma and Aeronomic Research (CSPAR), the Uni- 
versity of Alabama in Huntsville (UAH). Our code nu- 
merically integrates Eq. The s patial descritization 
employs a pseudosp e ctral algorithm ((Gottlieb et al 19771 : 
IShaikh fc Zank 20061 : IShaikh fc Zank 20071 ') based on a 
Fourier harmonic expansion of the bases for physical vari- 
ables (i.e. the magnetic field, velocity), whereas the tempo- 
ral integration uses a Runge Kutta (RK) 4th order method. 
The boundary conditions are periodic along the x,y and z 
directions in the local rectangular region of the solar wind 
plasma. 
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Figure 1. 3D simulation of whistler wave turbulence in the 
kde < 1 regime exhibits a Kolmogorov-like inertial range power 
spectrum close to fc~^/'^. The simulation parameters are: Box size 
is Lx X Ly X = 27r X 27r X 2it, numerical resolution is Nx x Ny X 
= 200x200x200, electron skin depth is de = 0.015, magnitude 
of constant magnetic field is Bq = 0.5. The characteristic large 
scales in this regime possess strong dispersion and wave activity 
that can be quantified from the turbulent equipartition between 
the velocity and magnetic fields. 



The turbulent fluctuations are initialized by using a 
uniform isotropic random spectral distribution of Fourier 
modes concentrated in a smaller band of lower wavenum- 
bers (fe < 0.1 kmax)- While spectral amplitudes of the fluc- 
tuations are random for each Fourier coefficient, it follows 
a certain initial spectral distribution proportional to k^", 
where a is an initial spectral index. The spectral distribu- 
tion set up in this manner initializes random scale turbulent 
fluctuations. We note that a constant magnetic field is in- 
cluded along the z direction (i.e. Bo = Bqz) to accommo- 
date the large scale (or the background solar wind) magnetic 
field. The size of the 3D computational domain is (27r)'^ with 
the spectral resolution 256'^. In this paper, we present the 
results of freely decaying whistler wave turbulence and fo- 
cus primarily on understanding the inertial range cascades 
in both the kd^ < 1 and kde > 1 regimes. In principle, 
turbulence can be driven. The driven whistler turbulence is 
nonetheless beyond the scope of this paper. 

Electron whistler fluid fluctuations, in the presence of 
a constant background magnetic field, evolve by virtue of 
nonlinear interactions in which larger eddies transfer their 
energy to smaller ones th rough a forward cascade. Accord- 
ing to l|Kolmogorov 194l] l. the cascades of spectral energy 



Figure 2. The small scales magnetic field fluctuations in the 
kde > 1 regime depicts a Kolmogorov-like fc~^/^ spectrum which 
is a characteristic of hydrodynamic fluid. The simulation param- 
eters are same as those used in Fig 1, except de = 0.15. Our 
simulations show that the small scale fluctuations evolve as non 
magnetized eddies of hydrodynamic fluid where whistler waves do 
not influence the energy cascades. 

occur purely amongst the neighboring Fourier modes (i.e. 
local interaction) until the energy in the smallest turbulent 
eddies is finally dissipated gradually due to the finite dis- 
sipation. This leads to a damping of small scale motions. 
By contrast, the large-scales and the inertial range turbu- 
lent fluctuations remain unaffected by direct dissipation of 
the smaller scales. Since there is no mechanism that drives 
turbulence at the larger scales in our model, the large-scale 
energy simply migrates towards the smaller scales by virtue 
of nonlinear cascades in the inertial range and is dissipated 
at the smallest turbulent length-scales. The spectral trans- 
fer of turbulent energy in the neighboring Fourier modes in 
whist ler wave turbulenc e follows a Kolmogorov phenomenol - 
ogy (jKolmoeorov 194ll : llroshnikov 19631 : iKraichnan 1965h 
that leads to Kolmogorov-like energy spectra. We find from 
our 3D simulations that whistler wave turbulence in the 
kde < 1 and kde > 1 regimes exhibits respectively k'"^^^ 
(see Fig 1) and k~^^^ (see Fig 2) spectra. The inertial 
range turbulent spectra obtained from our 3D simulations 
are further consisteri t with 2D work (Biskamp et al 19961 : 
iDastgeer et al. 2000al : lOastgeer et al. 2000bh . Interestingly, 
it is evident from the whistler wave dispersion relation that 
the wave effects dominate in the large scale, i.e. kde < 1, 
regime where the inertial range turbulent spectrum depictes 
a Kolmogorov-like k~'^^^ spectrum. On the other hand, tur- 
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bulent fluctuations in the smaUer scale (fede > 1) regime 
behave like non magnetic eddies of hydrodynamic fluid and 
yield a k~^^^ spectrum. The wave effect is weak, or neg- 
ligibly small, in the latter. Hence the nonlinear cascades 
are determined essentially by the hydrodynamic like in- 
teractions. The observed whistler wave turbulence spec- 
tra in the kd^ < 1 and kde > 1 regimes (Figs 1 & 
2) can be followe d from the K olmogorov-Uk e argu ments 
IIKolmogorov 19411 : llroshnikov "963 : Kraichnan 19651 ) that 
describe the inertial range spectral cascades. We elaborate 
on these arguments to explain our simulation results of Fig. 
(1) & (2) in the following section. 



4 ENERGY SPECTRA IN WHISTLER WAVE 
TURBULENCE 

The exact spectral indices corresponding to the whistler 
wave turbulent spectra, described by the ideal electron mag- 
netohydrodynamic invariant, can be un derstood from the 
Kolmogorov's dimensional argu ments (|Kolmogorov 194ll : 
llroshnikov 19631 : iKraichnan 1963 ). The electron skin depth, 
de, in EMHD turbulence intrinsically divides the entire 
Fourier spectrum into regions for which length scales are 
either larger or smaller than de. We derive the spectral in- 
dices for both regions of the turbulent spectrum. 

In the underlying whistler wave model of magnetized 
plasma turbulence, the inertial range eddy velocity is char- 
acterized typically by «e ~ V X B. Thus the typical velocity 
of the magnetic fleld eddy Bi with a scale size £ can be 
represented by Ve — Bi/£. The eddy turn-over time is then 
given by 

£ £^ 
r ~ — ~ — . 

Ve Be 

This is the time scale that predominantly leads to the non- 
linear spectral transfer of energy in fully developed whistler 
wave turbulence. While the inertial range nonlinear cascades 
are determined essentially by the eddy turn over or spectral 
transfer time scale, it is not clear whether the characteristic 
length scales bigger than de, where whistler wave propaga- 
tion dominate, are influenced by whistler interaction time 
scales. We will comment on this issue in the following along- 
with the inertial range spectra in both the kde < 1 and 
kde > 1 regimes. 

4.1 kde < 1: Whistler wave regime 

In the regime where characteristic length scales are bigger 
than the electron skin depth {kde < 1), the inertial range 
whistler turbulent energy is dominated by the large scale 
fluctuations. The total energy corresponding to the turbu- 
lent fluctuations in this regime is then given as, 

E ~ |B1^ ~ BI ~ Vei'^. 

The Be represent magnetic fleld associated with the mag- 
netic field eddy of length £. The second similarity follows 
from the assumption of an equipartition of energy in the 
magnetic and velocity fleld components of whistler waves. 
The process of equipartition origintes from the correla- 
tion between the velocity and magnetic fleld fluctuations 



Figure 3. Turbulent equipartion between the velocity and mag- 
netic fields is observed in our 3D simulations. The equipartition is 
measured for the entire turbulent spectrum at each time step by 
the relation Eequi — X/fc(l^'=l'^ — k'^\Bk\'^). When the characteris- 
tic turbulent modes evolve towards equipartition, the relationship 
|ve(k, i)P ~ A;^[B(k, i)p is obeyed. Consequently, Eequi 10^^. 
This number is small enough to establish a nearly perfect equipar- 
tition between the velocity and magnetic field associated with the 
whistler waves. 



We ~ k X B, where k = kxX + ky-y -I- kzZ is a three dimen- 
sional wave vector. The latter is further consistent with the 
electron flow speed, i.e. Eq (6) in combination with the wave 
perturbed magnetic fleld Eq. (11), that is used to derive the 
dynamical equation of whistler wave turbulence, i.e. Eq. (|8]). 
This velocity-magnetic fleld correlation essentially produces 
the velocity field fluctuations that are normal to the mag- 
netic fleld in a whistler wave packet. Consequently, the en- 
ergy associated with the velocity and magnetic field for each 
characteristic turbulent mode evolves toward a relationship 
that satisfies v1 ~ k^B^ . To quantify our arguments, we fol- 
low the evolution of turbulent equipartion in our simulations 
by computing the following quantity, 

Eequ.{t) ~ Y^{\Ve{k,t)\' - k^\B{k,t)f), (13) 
k 

which should be close to zero for the inertial modes that 
exhibit nearly perfect equipartition. The summation in Eq. 
(|13[l is carried over all the modes (i.e. fc's) that constitute 
the inertial range spectrum. Our simulation results, follow- 
ing the evolution of Eq. (|13|l . are shown in Fig 3. We find 
from our simulations that the inertial range turbulent fluctu- 
ations closely follow the equipartition that leads to a strong 
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wave activity in the fcde < 1 regime. Despite the pres- 
ence of the dispersive whistler waves in this regime, the 
inertial range spectrum continues to follow a Kolmogorov- 
like fc~^/^ spectrum. The equipartition in whistler wave 
turbulence h as also been reported in our two dimen- 
sional work jPastgeer et al. 2000al : IPasteeer et al. 2000bl : 
IShaikh fc Zank 20051 ). Interestingly, our 3D simulations, de- 
scribing the equipartition between the velocity and magnetic 
field fluctuations, are consistent with the 2D counterpart. It 
thus appears that the turbulent equipartition is a robust fea- 
ture of whistler waves that is preserved in both 2D and 3D 
nonlinear mode coupling interactions. The spectral cascades 
of inertial range turbulent energy is nonetheless determined 
by the energy cascade per unit nonlinear time as follows, 

~ £ ~ M 

On cissuming that the spectral energy cascade is loc al in the 
wavenumber sp ace (|Kolmogorov 1941; Kolmogo rov 194ll : 
llroshnikov 19631 ) . the energy spectrum per unit mode yields 
Ek ~ e" fc'' . On substituting the energy and energy dissipa- 
tion rates and equating the powers of Bg and £, we obtain 
Q = 2/3 and /3 = —7/3. This, in the fcde < 1 regime, leads 
to the following expression for the energy spectrum 

Ek ^ e^^''k-'^\ 

It can be noted from the dispersion relation, Eq. (|10|l . that 
the group velocity of whistler waves in the fcdg < 1 regime 
is ui/ky ~ duo/dky ~ k and the dispersion is uik ~ ujokyk. 
Both the quantities are proportional to the characteristic 
wavenumber k. It is evident from these relations that the 
group velocity and dispersion of whistlers are predominant 
at the smaller length scales in the kde < 1 regime. Cor- 
respondingly, the scaling law for energy cascades is modi- 
fied by the short scale spectrum of the whistler waves in 
the kde < 1 regime. To determine the effect of small scale 
whistler waves on the spectral transfer, we compute the en- 
ergy transfer rates in the kde < 1 regime as follows. 

The kde < 1 regime comprises the dispersive whistler 
waves whose interaction time can be estimated from r„ ~ 
£/vg, where Vg is the group velocity for the whistler modes. 
The group velocity of whistler waves in the kde < 1 regime 
is Vg ~ dio/dky ~ k ~ £~^. The interaction time between 
two (or more) whistler wave packets thus yields t^i ^ . 
The nonlinear energy cascade rates computed as above, i.e. 
e ~ E/t, will be modifled by the whistler interaction time 
as 

P- ■ 

Here is the whistler modified energy transfer rates. 
On using the Kolmogorov phenomenology that the spec- 
tral transfer is local and depends only on the en- 
ergy dissipation rates and mo des ([Kolmogorov 194ll : 
IKolmogorov 194ll : llroshnikov 19631 ). the energy spectrum 
can be given by Ek ~ s^k^. Upon substituting the en- 
ergy dissipation rates, we estimate the spectral energy as 
Ek ~ e^^^k~^. The change in the inertial range spec- 
tral slope due to the whistler w a ves is referred to as 
whistler effect (iBiskamp et al 19961: Dastgeer et al. 2000al : 



E 



iDastgeer et al. 2000bl : IShaikh fc Zank 20051 ). By introduc- 
ing the whistler time scale in deriving the energy cascade 
rates Ew, it is noteworthy that the spectrum in the kde < 1 



regime is modified by the presence of whistler waves and 
one might infer that the whistler waves modify the iner- 
tial range spectrum from k'"^^^ to a more flatter one, i.e. 
k^^. Although the difference between the two spectra is 
small enough to be noticeable in the 3D simulations (gen- 
erally because of poor spectral resolutions), the flattening 
of the spectrum is not observed in our simulations that 
persistently show that the whistler wave spectrum is close 
to fc"^''^. While the spectral resolution in our three di- 
mensional simulations is not adequate enough to resolve 
the two distinct sp ectra, very high resol ution simulations 
(upto 5120^) in 2D (|Shaikh fc Zank 20051 ) suggest that the 
volume integrated energy spectra are not affected by the 
presence of the whistler waves and the inertial range tur- 
bulent fluctuations continue to exhibit a Kolmogorov-like 



k 
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spe ctrum. The whistler effect in those simulations 



(Shaikh fc Zank 20051 ) is observed to be influential only at 
the local region in the inertial range turbulent spectrum. 
This r esult is further cons istent with that of MHD turbu- 
lence (jShebalin et al. 19831 ) where anisotropy in the spec- 
tral space mediated by the Alfven waves (i.e. the Alfven ef- 
fect) is explained by virtue of local Fourier mode, while the 
volu me integrated MH D spectrum exhibits a Kolmogorov- 
like ({K olmogo rov 194ll) k~^^^ pow er law. The controversy 
(llrosh nikov I9U iKraichnan 19651 ) with regard to the k 
or k~^^^ MHD spectrum is not a subject of this paper and 
we will not discuss it any furt her. The reader can however 
refer to the book by Biskamp (|Biskamp 20031 ) and the ref- 
erence therein. 



4.2 kde > 1: Hydrodynamic-Iike regime 

The regime kde > 1 in whistler wave turbulence corresponds 
essentially to a hydrodynamic regime because the EMHD 
equation, Eq. (|8)l, in this regime reduces to the Navier Stokes 
equation that describes the dynamics of non magnetized hy- 
drodynamic flows. The energy spectrum in this regime is 
dominated by the shorter length-scale turbulent eddies that 
give rise to the characteristic spectrum of an incompressible 
hydrodynamic fluid. The group velocity of whistlers, in this 
regime, is small and hence it is expected that the effect of 
whistlers will not be present. For kde > 1, the first term in 
the energy can be neglected and thus 



E ■ 



k'^B^ 



where we have used fc ~ l/£. The energy cascade rates per 
unit nonlinear transfer time, e ~ -E/r, in the regime kde > 1 
lead to e ~ Bf/£'^. On using the Kolmogorov's phenomenol- 
ogy of local spectral cascade, the energy spectrum of whistler 
turbulence in the kde > 1 regime can obtained as 



Ek 



_2/3 



fc" 



5/3 



in agreement with our simulations (see Fig 2). This spec- 
trum is identical to that of energy in three dimensional in- 
compressible Navier-Stokes turbulence and further confirms 
the hydrodynamic nature of the whistler wave turbulence 
for the small scale fiuctuations in kde > 1 regime. While the 
longer scales {kde < 1 modes) possess stronger tendency of 
behaving like whistlers, the shorter scales {kde > 1 modes) 
act like unmagnetized hydrodynamical eddies where wave 
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effects are considerably weaker. Hence whistler wave turbu- 
lence in this regime exhibits the energy spectrum that is 
essentially identical to that of hydrodynamic fluid. 



5 SUMMARY 

Three dimensional sinmlations of turbulent cascades in solar 
wind plasma arc carried out to quantify the role of whistler 
waves corresponding to the incrtial range fluctuations that 
possess characteristic frequency bigger than the ion gyro 
frequency (w > ajd) and length scales smaller than the 
ion gyro radius {kpi > 1). In this regime, the solar wind 
plasma fluctuations comprise of unmagnetized ions, hence 
the entire dynamics is governed by the electron fluid mo- 
tions. The rotational magnetic field fluctuations in the pres- 
ence of a background magnetic field lead to propagation of 
dispersive whistler waves in which the wave magnetic and 
velocity fields are strongly correlated through the equiparti- 
tion (Ug ~ k^B^). The latter is employed in our simulations 
to quantify the role of whistler waves that arc ubiquitously 
present in the inertial range in the high frequency {uj > cud) 
solar wind plasma. Interestingly we find that despite strong 
wave activity in the inertial range, whistler waves do not in- 
fluence the inertial range turbulent spectra. Consequently, 
the turbulent fluctuations in the inertial range are described 
by Kolmogorov-like phenomenology. Thus consistent with 
the Kolmogorov-likc dimensional argument, we find that 
turbulent spectra in the kde < 1 and kde > 1 regimes are 
described respectively by fc"'^^^ and k~^^^ . Our results are 
important particularly in understanding turbulent cascade 
corresponding to the high frequency (w > uid) solar wind 
plasma where characteristic fluctuations are comparable to 
the electron inertial skin depth. 

The support of NASA(NNG-05GH38) and NSF (ATM- 
0317509) grants is acknowledged. 
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